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An extensive photometric study of the Blazhko RR Lyrae 
star RZ Lyr* 
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ABSTRACT 

The analysis of recent, extended multicolour CCD and archive photoelectric, pho- 
tographic and visual observations has revealed several important properties of RZ Lyr, 
an RRab-type variable exhibiting large-amplitude Blazhko modulation. On the time- 
base of ~110 yr, a strict anticorrelation between the pulsation and modulation period 
changes is established. The light curve of RZ Lyr shows a remarkable bump on the 
descending branch in the small-amplitude phase of the modulation, similarly to the 
light curves of bump Cepheids. We speculate that the stellar structure temporally suits 
a 4:1 resonance between the periods of the fundamental and one of the higher-order 
radial modes in this modulation phase. The light-curve variation of RZ Lyr can be 
correctly fitted with a two-modulation-component solution; the 121 d period of the 
main modulation is nearly but not exactly four times longer than the period of the 
secondary modulation component. Using the inverse photometric method, the varia- 
tions in the pulsation-averaged values of the physical parameters in different phases 
of both modulation components arc determined. 

Key words: stars: horizontal branch - stars: oscillations - stars: variables: RR Lyr 
- stars: individual: RZ Lyr - techniques: photometric - methods: data analysis 



1 INTRODUCTION 

The light-curve modulation phenomenon of RR Lyrae stars, 
the so-called Blazhko effect, is still one of the unsolved prob- 
lems of stellar pulsation theory. Des pite the recen t groun d- 



2009bl. 2011; Sodor et al. 


12007. 2011; Chadid et all |2010|; 


GusKenberger et al. 1 201] 


; Kolenberg et al.l 2011T) and the- 


oretical (IStothersI 


20061. l20ld. 2011; Smolec et all l201ll; 


Szabo et al. 20ld; 


Kollath, Molnar & Szabd 2011) efforts 



aiming at disclosing the root cause of the light-curve 
changes, no convincing explanation for the phenomenon has 
been given yet. 

We have learned a lot about the complexity of the fre- 
quency spectra of the light variation, the occurrence rate of 
the modulation, the long-term changes of the modulation 



* Based on observations collected primarily with the automatic 
60-cm telescope of Konkoly Observatory, Budapest, Svabhegy 



properties and the changes in the mean global physical pa- 
rameters during the modulation cycle in recent years, but 
none of the proposed explanations is capable to elucidate all 
the observed features of these peculiar variable stars. There- 
fore, we think that further detailed studies may still help in 
understanding the Blazhko phenomenon. 

An ideal target for studying the Blazhko effect is RZ 
Lyr (< V >= 11.7 mag, P=0.511 d, a 2 ooo = 18 h 43 m 37?9, 
^2000 = +32° 47' 53'.' 9). Its photometric history dates back 



Williams! ( 1903 


) , the star has been the target of a large num- 


ber of visual ( 


Batvrevl 1951, 19521; Tsessevich 1953, 


1958; 


Klcpikova 1958: Bclik 1969: Romanov 1969; Bogdanov 


1972; 


Zverev & Makarcnkd 1 1979D and some photoelectric studies 



Il973l ; iButler et al.lll982f u" In addition, unprocessed, archive 
photographic and photoelectric observations of RZ Lyr are 
also available at the Konkoly Observatory. 
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Figure 1. Identification of the comparison stars and the new 
variables in the field of view of RZ Lyr. 

The large-amplitude, relatively long-period (~116 d) 
Bl azhko modu l ation of RZ Lyr was already established 
by iTsessevichl (|l955t ). Based on the quasi-monotonic pe- 
riod de crease of the star (a = P -1 dP/dt = —0.367 
Myr" 1 , |Le Borgne et al.ll2007al ), it could be supposed that 
the light-curve modulation of RZ Lyr does not show an 
erratic/complex behaviour either. Therefore, an extended 
CCD observation series of RZ Lyr was initiated in 2010 in 
order to follow the connection between the pulsation and 
modulation period changes of the star, and to disclose the 
changes in the mean values of the physical parameters dur- 
ing its Blazhko modulation. 



2 OBSERVATIONS OF RZ LYRAE 

The analysis of the pulsation and modulation properties of 
RZ Lyr is based mostly on a recent, extended CCD data 
set and on previously unpublished, archive photographic 
and photoelectric observations obtained with the telescopes 
of the Konkoly Observatory. In addition, all the published 
photometric information available on the light variation of 
RZ Lyr have been collected and are utilized. 



Table 1. Johnson-Cousins magnitudes of the photoelectric and 
CCD comparison stars 



Comp.star/2MASS ID B V I c Ref. 

a 18432371+3253475 11.546(19) 11.055(11) 10.454(29) AAVSO VSP a 

b 18434732+3256105 11.440(20) 10.814(10) 10.120(21) AAVSO VSP a 

c 18435761+3249546 12.846(16) 11.977(13) 11.009(18) AAVSO VSP a 

d 18443148+3257263 9.966 9.800 Sturch (19661 

a http://www.aavso.org-/vsp/ 

with the 60-cm telescope of the Konkoly Observatory in 
1958, 1959 with RCA 1P21, and in 1968, 1969 and 1972 
with EMI 9052 B photomultiplier tubes. About 700 obser- 
vations were obtained in B and V bands. Three different 
comparison stars were used during the observing runs. Stan- 
dard magnitudes of the comparison stars are taken from the 
comparison sequences of the AAVS O Variable Star Plottefl 
(5671dpv) and from ISturchl (|l966t l. The BVIc magnitudes 
of the comparison stars are summarized in Table [T] 

CCD observations of RZ Lyr were obtained on 109 
nights between 2010 April and 2011 September with the 
automatized 60-cm telescope of the Konkoly Observatory at 
Svabhegy, Budapest, equipped with a Wright Instruments 
750 x 1100 pixel CCD camera and BVIc filters. About 
2300 exposures were taken in each band. Relative magni- 
tudes of RZ Lyr were measured against the comparison star 
'a'. The photoelectric and CCD observations were trans- 
formed to Johnson-Cousins BVIc magnitudes by standard 
procedures. 

Originally, another blue star (2MASS 18431843 
+3248587) was selected as comparison for the CCD data, 
which, however, proved to be variable (VI). Without a sys- 
tematic variable search, another new variable (V2: 2MASS 
18435377 +3247459) was also identified in the field. VI and 
V2 were measured relative to the comparison stars 'a' and 
'c', respectively. The light curves of the new variables are 
given in the Appendix. The comparison stars and the new 
variables (V1,V2) are identified in Fig.fl] The Konkoly pho- 
tometric observations and maximum times of RZ Lyr are 
available as Supporting Information with the electronic ver- 
sion of this paper. Tables [2] and [3] give samples regarding 
their form and content. 



2.1 Konkoly observations 

A total number of about 1400 photographic exposures were 
taken with the 16-cm astrograph of the Konkoly Observa- 
tory at Budapest between 1950 and 1954, on 130 plates. The 
plates were digitalized using an EPSON Perfection V750 
flatbed transparency scanner. The photographic densities of 
RZ Lyr and the comparison stars were determined using 
standard IRAF0 aperture photometry packages. The magni- 
tudes of RZ Lyr w ere evaluated using Tycho-2 B magnitudes 
(|H0g et alj|200d ) of the surrounding stars. 

Photoelectric observations of RZ Lyr were collected 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which arc operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



2.2 Other photometric data 

In the GEOS3 (|Le Borgne et alj 12007*3 ) data base, all the 
published maximum times have been collected and many 
other maximum observations are also given. We have in- 
spected all the CEOS data by checking the given maximum 



available ( 


Belikl 


1969; Fitch, Wisniewski & Johnson 


1966; 


KleDikoval 


19581; 


Migach 


| 19691; |Romanov 19691; Tsessevichl 


19531. 1958 


; Tsessevich 1 


Il969l; IZverev & Makarenkol Il979h. 



Additionally, multiple entries are corrected and outlying 
points are removed. 

Besides maximum timings, light-curve data and 
maximum-brightness values are also used to detect changes 



2 http://www.aavso.org/vsp/ 

3 http://rr-lyr.ast.obs-mip.fr/ 
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Table 2. Konkoly photometric observations of RZ Lyr. The full 
table is available as Supporting Information with the online ver- 
sion of this paper. 



P,=0.51 1 247 E=2441 1 83.426 



HJD 

2400000+ 


mag 


Detector 


Filter 


55294.64087 


12.369 


CCD 


B 


55294.64436 


12.339 


CCD 


B 


55294.64238 


11.942 


CCD 


V 


55294.64341 


11.325 


CCD 


1 


36413.4595 


12.656 


pe 


B 


36413.4643 


12.145 


pe 


\ ' 


33369.501 


12.29 


Pg 





Table 3. Times and brightnesses of maxima redetermined from 
published visual and photoelectric light curves and determined 
from the Konkoly observations. The full table is available as Sup- 
porting Information with the online version of this paper. 



Times of maxima 


pg/B 


vis/V 


Ic 


references 


HJD-2400000 


mag 


mag 


mag 




26074.5430 




10.67 




1 


26209.4900 




11.09 




1 


33369.541 


11.04 






3 


33372.610 


11.06 






3 


36413.513 


10.98 


10.81 




3 


55346.356 


10.983 


10.842 


10.617 


3 


* reference numbers 


are the 


same as in 


Table H] 





in the pulsation and modulation properties. In Ta- 
ble [3] maximum times and brightness values deter- 
mined and/or rede termined fr o m the original lig h t-curv e 



data published in rTsessevichl (ll953Tl; iKlepikoval jl958h 
iFitch. Wisniewski fc Johnson! (1 19661 ): IZverev fc Makarenkol 
i|l979T ). which were utilized for determining the modulation 
periods (see Table [4} are also given. 

No good quality light curve of RZ Lyr was published 
between 1972 and 2010. However, thanks to the efforts 
of mainly the GEO 93 and BAV0 observers the maximum 
times available in the GEOS data base cover this pe- 



riod (lAgerer et al l l200l|; lAge 



2003: 



Huebs chcr 
20091. l20ld: 



rer fc Huebscherl Il99€ 
20051. 120111: iHuebscher et all 2005 _ 
Le Borgne et all 12005) [2006allbL l2007rJ lcL 



200i 

2008allbll3 . l2009al lbl. l20ld . l201ll : ISamolvlj[2010l . l201ll ) . 



200C 

2ooe 
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4 http://geos.webs.upv.es/ 

5 http://www.bav-astro.de/ 



Figure 2. Top panel: O — C values of the maximum-timing obser- 
vations of RZ Lyr. The O — C data reflect a quasi-monotonic pe- 
riod decrease during the ~110 yr of the observations; a parabolic 
fit to the O — C data is also shown. Different segments of the data 
are fitted with lines for a direct determination of the temporal 
period values. Bottom panel: pulsation and modulation periods 
calculated for different segments of the data. The temporal pul- 
sation periods that correspond to the O — C data segments are 
indicated by horizontal black lines, gray symbols show the pul- 
sation periods determined from the light-curve data. Modulation 
periods determined from maximum brightness, phase and light- 
curve analyses are plotted by filled and empty squares and filled 
circles, respectively. The linear pulsation and modulation period 
variations shown by straight lines correspond to the parabolic 
O — C fit and the linear fit of the Blazhko period values. 



3 LONG-TERM CHANGES IN THE 

PULSATION AND THE BLAZHKO PERIODS 

To follow the pulsation-period changes of RZ Lyrae, the 
O — C diagram is constructed using the maximum tim- 
ings discussed in Section T2.2I complemented with the maxi- 
mum times obtained from the Konkoly data (Table [3]). Ac- 
cording to the O — C variation shown in Fig. the pul- 
sation period has been decreasing more or less continu- 
ously during the ~110 yr covered by the observations. A 
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Table 4. Pulsation and Blazhko periods of RZ Lyr at different epochs 



Time-interval 


Po 


Error 


Diff* 


Pm 


Error 


Nmax 


Method** 


Data (ref)*** 


JD-2400000 


d 


io- 


7 d 


d 


d 








26074 - 32473 


0.5112501 


1 





116.1 


0.1 


88 


a 


visual (1), (2) 


33369 - 35041 


0.5112464 


4 


-11 


117.1 


0.2 


53 


a 


photographic (3) 


36413 - 42955 


0.5112419 


2 


-34 


118.04 


0.04 


34 


a 


photoelectric V (3), (4), (5), (6), (7) 


52416 - 54014 


0.5112462 


10 


72 


120.1 


0.5 


48 


1) 


maximum timings (8) 


54205 - 55826 


0.5112343 


7 


-39 


120.8 


0.2 


85 


b 


maximum timings (3), (8) 


55294 - 55826 


0.5112300 


1 


-80 


121.36 


0.04 




c 


CCD (3) 



* difference from the period value predicted from a parabolic O — C fit 

** a: pulsation and modulation periods from fitting light curve and maximum brightness data, respectively 
b: pulsation and modulation periods from fitting maximum timing data 

c: p ulsation and modu latio n periods from complete ligh t curve analysis 

*** (1 ) ITsessevichI <1953ft : (2) IZverev fc Makarenkol jl979j h (3) this paper; (4) ISturchl jl966t k f5l lFitch. Wisniewski fc Johnson! (l966t ): 
(6) IButler et al.l dl982T ); f 7) IRomanovl 1119731k (8) GEOS maximum timings 



parabolic fit to the O — C data yields an a — P 1 dP/dt = 
-8.85 x 1(T 10 d _1 (-0.323 Myr -1 ) period-decrease rate. 

However, the temporal periods defined by linear fits to 
the different segments of the O — C data indicate that the 
period change has not been strictly linear. Although differ- 
ences between the actually observed and predicted period 
values may arise from an incomplete averaging of the phase 
modulation of the Blazhko effect in the case of scarce data, 
this is not the case for the latest epochs. The maximum 
times available for the time intervals JD 2 452 416-2 454 014 
and JD 2 454 205-2 455 826 give complete phase coverage of 
the 120-d modulation cycle. The high accuracy of the max- 
imum timings collected by the GEOS team makes it also 
possible to detect the ~0.02 d phase modulation in these 
data sets. Therefore, besides the temporal pulsation peri- 
ods, the Blazhko periods are also derived for these epochs 
from maximum-timing data. 

The archive light curves and maximum-brightness val- 
ues enable us to determine both the pulsation and the 
Blazhko periods with high accuracy in a homogeneous man- 
ner for three epochs. In order to get the most reliable results, 
only the best-quality, most accurate and homogeneous light 
curves and maximum-brightness data obtained at the differ- 
ent time-intervals are utilized for this purpose. The pulsation 
periods are determined by fitting 6th-12th-order Fourier se- 
ries to the light-curve data, without taking the modulation 
components into account. The incomplete phase coverage 
of these data does not make a complete analysis possible. 
The Blazhko periods are obtained from sinusoidal fits to the 
maximum-brightness values. 

The first data set, which covers ~1 8 yr, com- 
prises the visual observat ions of ITsessevichI (|l953l ) and 
IZverev fc Makarenkol (| 19791 '). There is no significant system- 
atic difference between the magnitude scales of the two au- 
thors. These are the only available photometric information 
on the light variation of RZ Lyr for this period. Although 
a lot of visual observations are also available for the sec- 
ond half of the 20th century, the inhomogeneous quality 
and ill-defined magnitude scales do not render their use in 
the analysis possible. The pulsation and modulation periods 
are determined from the Konkoly photographic data for the 
1950-1954 interval. The third data set covering the 1957- 
1973 period includes the photoelectric observations pub- 
lished in this paper complemented with the photometric 



data oflSturchl Jl966l).lFitch. Wisniewski fc Johnson! (|l966h 
and IButler et al.l (|l982T ). Th e timings and m agnitudes of 
maximum light published by IRomanovl l|l973h are also uti- 
lized to derive the Blazhko period at this epoch. 

The recent CCD observations provide accurate pulsa- 
tion and modulation periods from detailed light-curve anal- 
ysis as discussed in Section 5. 

Table [4] summarizes the pulsation and modulation peri- 
ods derived for the different data sets using different meth- 
ods. 

In the bottom panel of Fig. [2J the pulsation and modu- 
lation period changes are shown according to the data given 
in Table [4] A strong anticorrelation between the pulsation 
and modulation period changes is detected; while the pulsa- 
tion period decreased by ~0. 000012 d the modulation period 
increased by ~5 d. A dPo/dP m = —2.4 x 10~ 6 period-change 
ratio is determined. While the modulation period values fit 
a linear period increase quite accurately, the pulsation pe- 
riod shows some erratic changes superimposed on the linear 
period decrease, especially at the latest epochs. 

In I Jurcsik et al.l (|201ll ). period-change data of Blazhko 
stars have been collected for all the variables with pul- 
sation and modulation periods known at least for two 
different epochs. The majority of these stars (10 out of 
13) show anticorrelated pulsation and modulation period 
changes while parallel changes in the periods are detected 
only in three cases (RR Gem, XZ Dra, V56/M5). All the 
variables that have enough data to determine both pul- 
sation and modulation periods for more than two epochs 
show, however, complex period-change behaviour. The cor- 
relation/anticorrelation between the pulsation and modula- 
tion periods manifests itself, in fact, as a tendency and not 
as a strict relation in these stars. Such a strict connection 
between the changes in the pulsation and modulation peri- 
ods as observed in RZ Lyr have not been detected in any 
Blazhko star previously. 



4 LIGHT-CURVE VARIATION DURING THE 
BLAZHKO CYCLE 

The CCD V light curves of RZ Lyr in ten different phases of 
the modulation are shown in Fig. [3j The amplitude variation 



Photometric study of the Blazhko star RZ Lyr 5 



10.6 F — r 




-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 
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Figure 3. Light-curve variation of RZ Lyr during the Blazhko cycle. The observations are plotted by gray dots; data belonging to a 
given Blazhko phase are highlighted in black; solid lines indicate the mean light-curve shape. The Blazhko phases (>!/) are given in the 
top left-side corners in the panels. 
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Figure 4. Light and colour curves in the small-amplitude phase 
of the modulation (*& = 0.5). A peculiar bump on the descend- 
ing branch is pronounced in this Blazhko phase. These light 
curves look very much like the light variation of bump Cephcids 
(eg. V Aql). 



is as large as 0.8 mag and the phases of the maximum vary 
within about 0.08 pulsation phase (~0.04 d). 

In Blazhko stars showing large amplitude modulations, 
the pulsation light curve is hevily distorted in the small- 
amplitude phase of the modulation, a s in the extreme case of 
V445 Lyr (|Guggenberger et al. |[201^ . in preparation), when 
the amplitude does not exceed 0.2 mag. In general, the am- 



plitude of the light curve is significantly smaller in these 
Blazhko phases than the amplitude of any non-modulated 
Blazhko star with the same pulsation period and metallicity 
value. 

The light curve of RZ Lyr is also the most anomalous 
at the small-amplitude phase of the modulation. A pro- 
nounced bump appears on the descending branch centred 
around pulsation phase 0.25-0.30, when the pulsation am- 
plitude is the smallest (see Fig. [4]|. No similar feature of 
any other Blazhko or stable-light-curve RR Lyrae star has 
been ever observed, the detected bumps in RR Lyrae light 
curv es do not appear earli er than phase 0.40 (see e.g. Table 
3 in lGillet fc Crowelll98Sl ). The double-peaked light curve, 
which shows a secondary maximum at the small-amplitude 
phase of the modulation of V455 Lyr (|Guggenberger et al. I 
|2012| . in preparation) may represent an extreme example of 
the bump phenomenon of RZ Lyr. 

The bump on the descending branch on the small ampli- 
tude light curve of RZ Lyr is reminiscent to the li ght curves 
of bump Cepheids jWisniewski fc Johnson! Il968l ). The pe- 
riods of the fundamental and the second-overtone radial 
modes are in 3:1 r esonance in these, double -mode, single- 
periodic pulsators jSmolec fc Moskalikll2010l ). resulting the 
bump in the light curve. The bump is also detected in the 
colour curves both in bump Cepheids and in RZ Lyr, in- 
dicating that the temperature-decrease is less steep during 
the bump phase. Therefore, we speculate, whether a simi- 
lar explanation for the appearence of the bump in RZ Lyr 
also holds. According to linear and non-linear nonadiabatic 
models (G. Kovacs and Z. Kollath private communications) 
the period of the fifth-overtone mode, which is, however, 
strongly damped in RR Lyrae stars, may be in 3:1 reso- 
nance with the fundamental-mode period in metal deficient 
variables if the luminosity is low enough (L ~ 35 — 501/0). If 
this were indeed the case, that is, a 3:1 resonance condition 
were temporally fulfilled in the low-amplitude phase of the 
modulation, then a peculiar behaviour of the amplitude and 
phase of the third-harmonic-order component of the pulsa- 
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Figure 5. Fourier amplitudes and amplitude ratios of the pulsa- 
tion light curves in different Blazhko phases. The amplitudes of 
the harmonic components show parallel changes with decreasing 
amplitudes in increasing harmonic orders. The normalized ampli- 
tude variations (top right-hand panel) indicate that the absolute 
strength of the amplitude variation decreases also with increas- 
ing harmonic orders. The A^/A-i amplitude ratios of the fourth 
and sixth harmonic components (filled symbols) are anomalously 
large at Blazhko phase 0.5 (bottom left-hand panel). This pecu- 
liar behaviour is even more evident if the A^^/A^ amplitude 
ratios are plotted (bottom right-hand panel). 



tion frequency would be expected to be detected when the 
bump is pronounced. 

Figs. [S] and [5] show the amplitude and phase variations 
of the six lowest-order pulsation components (/o,2/o, ..6/0) 
in ten different phases of the modulation. An anomalous 
behaviour of the Fourier components at the low-amplitude 
phases of the modulation (\& ~ 0.5) is indeed detected, how- 
ever, it is the most prominent in the amplitudes of the fourth 
and the sixth-order components (shown by filled symbols 
in the figures). Each amplitude ratio has a minimum value 
at ^ = 0.4 but with different depths. The minima of the 
variations of A2/A1, A3/A1 and A5/A1 are deeper than the 
minima of A4/A1 and Aq/A\. Moreover, A4/A1 and Aq/Ai 
have anomalously large values at * = 0.5 (bottom left-hand 
panel in Fig. [SJ. Consequently, the relative amplitudes of 
the consecutive orders A4/A3 and Aq/A$ have maxima at 
* = 0.4-0.5, while A2/A1, A 3 /A 2 and A5/A4 show minima 
in this phase (bottom right-hand panel in Fig. [SJ). 

As for the phases of the harmonic components of the 
pulsation frequency, $21 hardly varies relative to $1, while 
the higher-order phase differences show ~1 rad amplitude 
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Figure 6. Fourier phase of /o and phase differences (^k^i) of the 
higher-order harmonic components of the pulsation light curve in 
different phases of the modulation. The phase variation of the /o 
and 2/o components ($1 , $2.1 ) is only about 0.2 rad, while the 
phase variations of the higher harmonic-order components are 
as large as 1 rad. The amplitudes of the phase variations of the 
higher-order components are somewhat different as shown in the 
right-hand panel, which plots the absolute values of the phase 
difference variations. The amplitude of the variations in $41 is 
the smallest. 



variations with a maximum at ^ = 0.4-0.5. The amplitudes 
of the higher-order phase-difference variations differ slightly, 
$41 exhibits the smallest-amplitude variation. 

Consequently, if resonance causes the appearence of the 
bump, then it is most probably a 4:1 and not a 3:1 reso- 
nance. Recently, a 9:2 (4.5:1) half integer resonance with 
higher-order (8th-10th) radial modes have been suggested 
to explain the period doubling behaviour of s ome Blazhko 
stars observed by the Kepler space t elescope ()Szab6 et al.l 
l2010l : iKollath. Molnar fc SzaboTl201ll ). These higher-order 
modes have also been shown to be less damped as the third- 
seventh order modes. Their periods can match the 4:1 inte- 
ger resonance condition in certain cases, as well. Therefore, 
we propose that a 4:1 resonance with a higher-order radial 
mode explains the bump-shape character of the light and 
colour curves of RZ Lyr during the small-amplitude phase 
of the modulation. The feasibility of such a scenario has to 
be checked, however, by detailed hydrodynamical modelling. 

Another, natural explanation for the bump would be a 
shock- wave propagation through the atmosphere, as shock 
waves are responsible for the appearence of the bumps and 
humps observable at around 0.7 and 0.9 pulsation phases 
(|Gillet fc Crowdll98Sl : lc"hadid, Gillet fc Fokinll2000l ). which 
are common features of both Blazhko and non-Blazhko RR 
Lyrae stars. However, no hydrodynamic study of the pulsa- 
tion of RR Lyrae stars points to the existence of any shock 
wave at ~0.3 pulsation phase, when the stellar radius is the 
largest. 



5 FREQUENCY ANALYSIS OF THE CCD 
OBSERVATIONS 

The freque ncy anal ys is is performed using the MUFRAN 
package (|Kollathl |l99Ch , the CLEAN Algorithm 
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Figure 7. The spectral window, the dirty and cleaned spectra of 
the data prcwhitcned for the pulsation and the main modulation 
components (fc/o, fc/o±/m) obtained with the CLEAN Algorithm 
are shown in the top, middle and bottom panels, respectively. Al- 
though the daily alias components dominate the window struc- 
ture, based on the cleaned spectrum the residuals can be identified 
as a secondary modulation, as the frequency components are lo- 
cated symmetrically at both sides of the main pulsation frequency 
components. 



ijRoberts. Lehar fc Dreheilll9871 ), the GNUPLOT0 utilities 
and a nonlinear fitting algorithm developed by A Sodor. To 
a first approximation, the frequency spectrum of the data 
can be described by frequency triplets corresponding to 
the ~121-d (/ m = 0.0082 d _1 ) modulation of the pulsation 
light curve. 

A detailed analysis of the residuals reveals several ad- 
ditional frequency components. Fig. [7] shows the spectral 
window and the dirty and cleaned Fourier amplitude spec- 
tra of the data prewhitened with the pulsation and the main 
modulation components (triplets) . According to the cleaned 
spectrum, the residual can be described as a secondary mod- 
ulation with components at around ±4/ m separations from 
the main pulsation components. 22 frequencies matching 
this secondary modulation sequence have amplitudes larger 
than the 3.5 a level of the residual spectrum in at least one 
band; among these, the ratios of the observed amplitudes 
over their formal errors are larger than 4.0 for 21 compo- 
nents in at least two bands. 

15 pulsation frequency components are detected, how- 
ever, 14/o and 15/o are only marginally significant with am- 
plitudes corresponding to 3.0, 2.0, 1.8 and 3.8, 1.9, 2.3 times 
the mean level of the residual spectra in the B, V and Ic 
bands, respectively. 25 and 22 modulation components be- 
longing to the ±/ m and ±4/ m sequences are identified; / m 
and 4/ m are also detected. A long-term drift (/' ~ 0.002 
d _1 ) remains in the residuals after prewhitening the data 
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Figure 8. Schematic plot of the CCD V light-curve solution. The 
detected frequencies are indicated by vertical lines, with length 
corresponding to their amplitudes, in the vicinity of the pulsation 
frequency and its harmonics (fc/o). Note the different magnitude 
scales of the plots. The arrows indicate that the amplitudes of 
/o,2/o,3/o and 5/n are out of the ranges. Several modulation com- 
ponents are present in each harmonic order, the most prominent 
features arc the triplets with / m = 0.0082d~ 1 (121 d) separation. 



with the above frequencies, which has, most probably, in- 
strumental instead of stellar origin. 

After these prewhitening steps, the residuals still show 
the largest-amplitude signals in the vicinity of the pulsa- 
tion frequencies. Signals appear in these spectra at around 
fc/u ± 2/ m , ±3/ m and +5/ m ; the amplitudes of 26 compo- 
nents belonging to these series are larger than the 3.5<r level 
of the residual spectrum in at least two bands, while in 6 
less certain cases the amplitudes are larger than the 3. 0<t 
level in one band. 
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Table 5. Independent frequencies of the CCD light-curve solution 




frequency 

Figure 9. Final residual spectra in the B, V and Iq bands and 
their smoothed, 3<r levels defined as 3 times of the mean level of 
the residual spectrum are shown. There is no clear evidence of 
any additional frequency component. 



Finally, 11 components at around ±8/ m separation from 
the pulsation components are detected with 3.5a level of 
the residual spectrum in at least two bands, while in 5 less 
certain cases the amplitudes are larger than the 3.0a level 
in one band. 

Fig. [8] delineates the result of the light-curve anal- 
ysis schematically. The detected frequencies around the 
pulsation components are drawn by their amplitudes ob- 
tained for the V light-curve solution. The amplitude dis- 
tribution of the modulation components shows a unique 
feature. The amplitudes of the frequencies at ±4/ m sep- 
aration are significantly larger than the amplitudes of 
the ±2/ m components, moreover, the frequencies at ±8/ m 
separations have similar amplitudes as the ±2/ m quintu- 
plet components. This is in high contrast with the re- 
sults of all the detailed analyses of Blazhko light curves 
where quintu plets and/or higher order multiplet compo- 
nents appear (IHurta et al.ll2008l; I Jurcsik et al.ll2008al.l2009al; 
Poretti et a l. 2010; Chad id et al.ll201(J ; lGuggenberger et al. I 
20111 ; iKolenberg et alj|201ll ). The amplitudes of the multi- 
plet components typically decrease with increasing multi- 
plet orders according to the analysis of all the well observed 
Blazhko stars' light curve. Another general feature of the 
modulation is that, if higher multiplet-order components are 
detected, then the entire 'spectrum' of the multiplet is more 
or less equally populated. These properties are valid espe- 
cially in the lower pulsation orders. 

The discrepant behaviour of the amplitudes of the mul- 
tiplets observed in RZ Lyr suggests that the ±4/ m and ±8/ m 
components belong to an independent modulation quintu- 
plet. The components of this secondary modulation are close 
to the 4th and 8th order multiplets of / m but do not exactly 
equal with these frequencies. The residual scatter of the 
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Figure 10. Maximum-brightness and maximum-phase variations 
according to the CCD V maximum observations given in Table[l0] 
First and fourth-order fits to the data are also drawn. The fourth- 
order components of the fits have well pronounced amplitudes 
indicating that, beside the 121-d primary modulation, a ~30-d 
periodicity in the light-curve variation is also present. 



light-curve solution corresponding to this scenario is smaller 
by ~10 per cent than the rms of the one- modulation solution 
as documented in Table [5] Adopting two independent mod- 
ulation components (/ m i and / m 2), the frequency solution of 
the V data yields a frequency value for the secondary mod- 
ulation that differs by more than 30<r from the exact ±4/ m 
value. 

It is also important to note that the amplitudes of the 
fc/o — /mi negative components of the main triplets are larger 
than the amplitudes of the fc/o + /mi positive components, 
while the secondary modulation shows an opposite character 
(fc/o + /m2 > fc/o — /m2). The dissimilarity of the / m i and 
the ~ 4/ m i = / m 2 modulations supports interpreting them 
as independent modulations. 

The maximum-brightness and the maximum-phase val- 
ues of the CCD data (Table [3} are phased with the 121-d 
modulation period in Fig. 1101 First and fourth-order har- 
monic fits to the data are drawn in the Figures. Unam- 
biguously, the amplitudes of the fourth-order components 
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Table 6. Frequencies, amplitudes (in mag) and phases (in rad) of the B, V and Iq light-curve solutions of RZ Lyr. The Ao values are 
the magnitude zero points relative to comparison star 'a'. The full table is available as Supporting Information with the online version 
of this paper. 



Identification 


/ [cd- 1 ] 


B 

Amplitude 


Phase 


V 

Amplitude 


Phase 


Ic 

Amplitude 


Phase 


A 




0.548 




0.658 




0.696 




l/o 

2/0 
3/0 


1.956067 
3.912133 
5.868200 


0.5110(4) 
0.2000(4) 
0.1382(4) 


1.521(1) 
5.305(2) 
2.884(3) 


0.3931(3) 
0.1582(3) 
0.1110(3) 


1.481(1) 
5.313(2) 
2.885(3) 


0.2395(3) 
0.0982(4) 
0.0720(4) 


1.328(1) 
5.273(4) 
2.879(5) 


/ml 
2/ml 

fm2 
I/O — /ml 
2/0 — /ml 


0.008243 
0.016486 
0.033330 
1.947824 
3.903890 


0.0042(4) 
0.0031(4) 
0.0045(4) 
0.0903(4) 
0.0554(4) 


1.9(1) 

3.3(1) 

4.5(1) 

3.496(4) 

0.978(8) 


0.0026(3) 
0.0026(3) 
0.0038(3) 
0.0704(3) 
0.0448(3) 


1.5(1) 

3.6(1) 

4.64(9) 

3.511(4) 

0.992(7) 


0.0012(4) 
0.0021(4) 
0.0042(4) 
0.0452(3) 
0.0284(4) 


0.5(3) 

3.7(2) 

4.6(1) 

3.547(8) 

1.02(1) 


1/0 + /ml 
2/0 + /ml 


1.964310 
3.920376 


0.0570(4) 
0.0209(4) 


0.273(7) 
4.34(2) 


0.0437(3) 
0.0171(3) 


0.295(7) 
4.36(2) 


0.0275(3) 
0.0105(4) 


0.30(1) 
4.34(3) 


I/O — /m2 
2/0 — /m2 


1.922737 
3.878803 


0.0097(4) 
0.0121(5) 


2.93(4) 
0.41(4) 


0.0059(3) 
0.0087(3) 


2.92(6) 
0.47(4) 


0.0018(4) 
0.0053(4) 


3.0(2) 
0.47(7) 


1/0 + /m2 
2/0 + /m2 


1.989397 
3.945463 


0.0255(5) 
0.0202(4) 


4.00(2) 
1.39(2) 


0.0190(3) 
0.0168(3) 


4.06(2) 
1.42(2) 


0.0118(4) 
0.0111(4) 


4.16(3) 
1.55(4) 



of both the maximum-brightness and maximum-phase data 
are significantly larger than the amplitudes of the second 
and third harmonic-order components. This result confirms 
that the modulation is not singly periodic, and that the pe- 
riod of the secondary modulation is very close to P m /4. 

We thus conc lude that RZ Lyr, similarly to CZ Lac 
(|S6dor et al.l l201lT ). exhibits two modulations with modu- 
lation frequency values being very close to the 1:4 ratio. 
Table [6] gives the frequencies of the light-curve solution and 
the amplitudes and phases of the frequencies for the B, V 
and Ic time series. The complete table is available from the 
electronic version of this paper. The light variation of RZ 
Lyr is described by three independent frequencies /o,/mi 
and / m 2- In this representation, the modulation frequencies 
at around ±3/ m and ±5/ m are fitted as combinations of the 
two modulations, i.e. k(—f m i + / m 2) and fc(/ m i + /ma). 

The rms scatters of the full light-curve solutions that 
indicate both the goodness of the fits and the accuracy of 
the observations are 0.012, 0.009 and 0.010 mag for the B, V, 
and Ic data, respectively. The residual spectra and their 3a 
levels, which are calculated as 3 times the smoothed mean 
values of the residual spectra, are shown in Fig. [9] for the 
three bands. 



6 PHYSICAL PARAMETERS AND THEIR 

VARIATIONS DURING THE MODULATION 

The variations of the pulsation-averaged atmospheric pa- 
rameters of RZ Lyr during the Blazhko cycle are derived us- 
ing the inverse photometric B aade-Wesselink method (IPM; 
ISodor. Jurcsik fc Sz cidl 2008). These parameters are the ef- 
fective temperature (T e ff), luminosity (L), radius (R) and 
effective surface gravity (logg). The IPM needs only mul- 



ticolour photometric time s eries and synthetic colou rs from 
static atmosphere models (|Castelli fc Kuruc3 120031 ) as in- 
put. We applied this me thod already for fiv e modulated 
RRab stars (MW Lyr - lJurcsik et all l2008bl; DM Cyg - 
lJurcsik et aLll2009al ; RR Gem, SS Cnc - IS6dorll2009l ; CZ Lac 
201 ll) succesfully. 



Sodor et al 



6.1 Constant parameters 

In the first step, the Blazhko-phase independent parameters 
of RZ Lyr have to be determined. These are the metallicity 
([Fe/H]), mass (SOT) , distance (d) and the interstellar red- 
dening [E(B-V)]. 

The metallicity of RZ Lyr was measu red by 
ISuntzeff. Kraft fc Kinmanl (| 19941) (AS = 9.89) and iLavdenl 
(|l994l ) ([Fe/H1= -2.13). On a combined, homogenized 



metallicity scale (|jurcsikl 1 19951 ; lJurcsik fc Kovacsl Il996l ). 
these observations correspond to —1.84 and —1.85 [Fe/H] 
values, respectively. Thus, we accept [Fe/H]= —1.85 for the 

metallicity of RZ Lyr. 

As it was shown in lJurcsik et al.l (|2009bl ). the well- 
defined mean pulsation light curves of Blazhko stars yield 
correc t metallicity value accor ding to the [Fe/H](P, $31) for- 
mula (|Jurcsik fc Kovacslll996l ). This is not the case, however, 
for RZ Lyr. The photometric metallicity defined by the mean 
V light curve is —1.44. The anomalous light-curve shape at 
the small-amplitude phase of the modulation (see details in 
Section [3} may account for the incorrect photometric metal- 
licity estimate. If data belonging to Blazhko phases 0.4-0.6 
(the light curves in these phases show an anomalous bump 
on the descending branch) are removed from the data set, 
then the photometric metallicity is —1.65. This is closer to 
the spectroscopic metallicity value, but it is still 0.2 dex 
more metal-rich than expected. We thus conclude that the 
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Table 7. The mean physical parameters and their uncertainties are derived from the BVIq mean pulsation light curves of RZ Lyr by 
the inverse photometric method (IPM) using static atmosphere models with [Fe/H] =—1.85 metallicity. The luminosity is selected to 
match the evolutionary possible values for a metal-deficient horizontal-branch star and the mean physical parameters are required to fit 
the pulsation period via the pulsation equation. 



<M d M v R L T eS 

<m e (pc) (mag) (R Q ) (L ) (K) 



0.66 ±0.04 1380 ±60 0.739 ±0.004 5.12 ±0.20 47.5 ± 2.5 6670 ± 150 



applicability of the metallicity formula for the light curves 
of Blazhko stars is questionable in some cases, even if the 
mean light curve is well-defined. 

The interstellar reddening towards the direction of 
RZ Lyr is 0.09 mag according to t he extinction maps of 
ISchlegel. Finkbeiner fc Davis! l|l99Sft . iBrancol (|l992l ) derived 
E(B — V) — 0.07 mag for RZ Lyr from the near-minimum- 
light B — V colour. The average of these two estimates, 
0.08 mag, is accepted for the reddening of RZ Lyr. The 
unreddened colours of RZ Lyr, which also depend on the 
uncertainties of the magnitudes of the comparison star (Ta- 
ble [TJ, are supposed to be accurate within about 0.03 mag. 

The next two constant parameters we have to determine 
are the mass (9Jt) and the distance (d). In principle, IPM is 
capable to drive these parameters by the application of the 
method on the mean pulsation light curve. 

However, if the fitting process is allowed to adjust these 
parameters unconstrained, then the resulting mass and lumi- 
nosity values are too high if compared to horizontal-branch 
(HB) evolutionary model results. Maybe, as a result of the 
bump-shaped light curve at the minimum-amplitude phase 
of the modulation, the mean light curve of RZ Lyr re- 
mains also somewhat anomalous. Therefore, the luminos- 
ity is constrained to match the oxyg en-enhanced HB evo- 
lutionary models of iDorma 3 l| 19921 ). The possible lumi- 
nosity range of the insta bility strip of the [Fe/H] = — 1.78 
HB models is 45-57 L jDormanl fl992l. fig. 4). A modi- 



fied pulsation equation (|S6dor. Jurcsik fc Szc idl 2008, eq. 2) 
is applied to link the mass, radius, temperature and pe- 
riod values during the fitting process of the IPM. These 
constraints yield 9)1 — 0.62 -0.83 DJIq for the possible mass 
range. The HB evolutionary models contradict, however, 
the solutions with masses above 0.7 OT©, as these model 
tracks do not enter the instability strip at the given metallic- 
ity. Therefore, L = 47.5 ± 2.5 L© and DJl = 0.66 ± 0.04 9Jt 
are accepted for the possible luminosity and mass ranges 
of RZ Lyr. These parameters result in a d = 1380 pc dis- 
tance estimate with ±60 pc uncertainty depending on the 
accepted values/uncertainties of the luminosity, the V light 
curve magnitude zero point and the interstellar absorption 
[Ay = 3.14E(B - V)]. 

Applying these constraints, the goodness of the IPM fits 
are only marginally worse than in the let-free case. 

The pulsation- and modulation-phase-averaged mean 
radius, temperature and absolute visual brightness (Mv) 
of RZ Lyr are derived applying the IPM for the BVIc 
mean pulsation light curves with fixed luminosity, mass and 
distance values corresponding to the above discussed con- 
straints. 

Similarly to the previous analyses, the IPM code is 
run with four different internal settings (for details see 



ISodor. Jurcsik fc Szeidll |2008l , table 1) to estimate the in- 
herent uncertainty of the method. 

The results for the mean physical parameters of RZ Lyr 
are summarized in Table [7] The error ranges are estimated 
by taking into account the inherent error of the method, 
the uncertainties of the input parameters and the possible 
ranges of the constrained parameters. 



6.2 Blazhko-phase dependent parameters 

Fixing the Blazhko-phase independent physical parameters 
to their values determined in Section 16.11 the IPM is run 
using the light curves of RZ Lyr in 10 different phases of 
the modulation. In this way, the radius, luminosity, tem- 
perature and surface-gravity changes during the pulsation 
in different phases of the modulation are determined. The 
pulsation-averaged mean values of the variations of these pa- 
rameters characterize the global mean changes in the stellar 
parameters during the modulation cycle. Similarly to the 
IPM results of the previous investigations of Blazhko stars, 
the actual values of the fixed or constrained physical parame- 
ters influence only the averages of the other physical param- 
eters, and affect their variations during the Blazhko cycle 
marginally. A small amplitude difference between the solu- 
tions is detected if the input parameters are varied within 
their possible ranges (Table[7]), but the relative sign of their 
variations is unaffected. 

As RZ Lyr has two independent modulation periods, 
and the shorter one is approximately, but not exactly, one 
quarter of the longer one, the mean physical parameter 
variations are investigated both on the original data sup- 
posing a single-periodic modulation, and on synthetic data, 
which separate the primary and secondary modulation com- 
ponents. The separation of the modulations is carried out 
by generating synthetic light curves from the light-curve 
solution given in Table [H] using pulsation frequencies and 
modulation terms belonging to only one of the modulations, 
namely ±/ m i, ±2/ m i and ±3/ m i in one case, and ±/ m 2 and 
±2/ m 2 in the other case. The minute coupling modulation 
terms are omitted, as they are not periodic with any of the 
modulations. Considering their rather small amplitudes, the 
omission of these terms has a negligible effect on the re- 
sults (a simil ar treatment was applied in the IPM analy- 
sis of CZ Lac (|S6dor et al.ll201lh . another multiple-periodic 
Blazhko star). 

The IPM results of the original time series are plotted 
against the phases of the 121. 4-d dominant modulation in 
Fig. 1111 As a result of the secondary modulation, which has 
a period approximately one quarter of the primary modu- 
lation period, the resulting curves have, in principle, wavy 
shapes but the 10 phase bins of the main modulation are 
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Figure 11. Variations of the mean observed (left-hand pan- 
els) and the derived physical parameters (right-hand panels) of 
RZ Lyr in different phases of the primary (~121 d) modula- 
tion. The secondary (~30 d) modulation has not been removed 
from the data, this modulation accounts partially for the scat- 
ter e.g. in the magnitude averaged V brightness and the radius 
and log (/ variations. Magnitude- and intensity-averaged bright- 
nesses and colours are denoted by angle and round brackets, re- 
spectively. From top to bottom, the left-hand panels show: am- 
plitude modulation - total pulsation amplitude; pulsation-phase 
or period modulation - variation of the phase of the / p pulsa- 
tion component (empty circles) and deviation of the instanta- 
neous pulsation period from the mean pulsation period (contin- 
uous line); pulsation-averaged V magnitude; pulsation-averaged 
B — V colour; pulsation-averaged V — Iq colour. From top to 
bottom, the right-hand panels show the pulsation averages of the 
following physical parameters: radius; surface gravity; absolute 
visual magnitude; luminosity; effective surface temperature. 



insufficient to resolve these features. Therefore, most of the 
phase-bin-to-phase-bin variations (these are the most promi- 
nent in the magnitude-averaged brightness, radius, and log g 
variations) are probably real features connected to the sec- 
ondary modulation and not the consequences of the inherent 
noise of the method and/or the inaccuracy of the input fight 
curves. 

In Fig. 111! the observed (feft-hand panefs) and de- 
rived (right-hand panels) variations of the pulsation-period- 
averaged parameters during the Bfazhko eyefe of the main 
modulation are plotted. The observed parameters (from top 
to bottom) are the fulf ampfitude, the $1 phase variation of 
the main pulsation frequency (/o), and the different averages 
of the fight and cofour curves according to magnitude and 
intensity-scafe representations. The sofid fine drawn in the 



phase-variation panef indicates the pufsation-period varia- 
tion corresponding to the observed phase variation of fo- 
The totaf range of the period variation is about 0.0006 d 
during the Blazhko cycle. No significant phase fag between 
the ampfitude and period variation of the puisation fight 
curve is detected during the main, 121-d modufation. 

The mean physical-parameter variations shown in the 
right-hand panels of Fig. [11] indicate an extremely farge (6.3 
per cent) luminosity variation and modest, 0.8 and 0.4 per 
cent radius and temperature changes, respectiveiy. The vari- 
ations in the mean radius, iuminosity and temperature are 
aif in fine with the ampfitude variation of RZ Lyr, that is, 
these parameters have maximum vaiues when the pufsation 
amplitude is the largest. 

The IPM sofutions of the synthetic data show the vari- 
ations of the physical parameters along the two modulation 
components separately (Fig. I12[) . The results for the main 
modulation component agree welf with the resuits obtained 
for the originai fight-curve data (Fig. Hip , while the physicai- 
parameter variations aiong the secondary modufation show 
a different behaviour. The ampfitude variation is the oppo- 
site of the phase variation in this modufation, that is, when 
the ampfitude is the iargest the period is the shortest. No de- 
tectabfe temperature variation and very smaif (0.2 per cent) 
luminosity change during the ~30-d modulation are derived. 
The 0.01 Rq radius variation, which is in anti-phase with the 
amplitude variation and is in phase with the period variation 
is the most prominent change connected to this modulation 
component. It is also worth to note that ampfitude variation 
dominates the modulation of the primary and phase modu- 
lation of the secondary modulations. The range of the period 
variation in the secondary modufation is iarger compared to 
the range of the period variation during the primary modu- 
fation. 



7 SUMMARY 

Although the number of Blazhko stars with observations 
alfowing detailed anaiyses is increasing, we can stiff fearn a 
fot about the modulation phenomenon from these data. The 
anaiysis of RZ Lyr disciosed some novefties of the Bfazhko 
modulation, too. 

• Based on multi-epoch observations, a strict anticor- 
reiation between the pulsation and modulation period 
changes is found with dP /dP m = —2.4 x 10~ 6 or 
(dP /fb)/(dP m /P m ) = -5.6 x KT 4 ratio. 

There is a tendency that the pulsation and modufation pe- 
riods of Bfazhko stars exhibit preferabiy opposite-direction 
changes, however, no strictfy monotonic and stabie period- 
change behaviour has been detected previousiy. The strong 
anticorrelation between the periods of RZ Lyr suggests that 
the modulation period cannot have an ad hoc value, its 
fength and its variation is defined by the steifar parameters 
and their variations, as in the case of the puisation period 
vaiue. 

• A pronounced bump on the descending branch of the 
fight curve in the smaif-ampfitude phase of the modulation 
is detected. Based on the similarity of the appearence of 
the bump to the light-curve shapes of bump-Cepheids we 
speculate that some resonance mechanism is responsible for 
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Figure 12. Variations of the physical parameters of RZ Lyr according to the phases of the primary (~121 d) and secondary (~30 d) 
modulations. The two modulations are separated using the appropriate frequencies of the light-curve solution given in Table [6] The 
panels are the same as described in the caption of Fig. 1111 



the bump phenomenon in this case as well. Without a de- 
tailed hydrodynamical modelling, which is, however, out of 
the scope of the present paper, this possibility can neither 
be proved nor rejected. 

No RR Lyrae star with stab le light curve is known that 
shows either period doubling (|Szab6 et alJl2010T ) or a bump 
on the descending branch of the light curve at around 0.3 
pulsation phase. Supposing that 4.5:1 and 4:1 half- integer 
and integer resonances are responsible for the intermittent 
occurrence of these phenomena in Blazhko stars, it can be 
concluded that the resonance conditions never come true in 
a normal RRab star. These are unique, temporal specialities 
of some Blazhko stars at certain phases of the modulation. 
If this is indeed the case, then it means that the structure 
of a Blazhko star may differ significantly (at least in some 
Blazhko phases) from that of a normal RR Lyr star. The 
detected changes of the mean physical parameters during 
the Blazhko cycle (as the 6 per cent change in the mean 
luminosity) may support such a scenario. 

The onset of period doubling is not strictly connected 
to any Blazhko phase in RR Lyr but the bump feature ap- 
pears only in the small-amplitude phase of the modulation 
in RZ Lyr. The modulation period of RR Lyr is ~ 40 d, sig- 
nificantly shorter than the 121-d Blazhko period of RZ Lyr. 
Moreover, the modulation itself is highly variable in RR Lyr 
while it is rather regular in RZ Lyr. Therefore, it is a plau- 
sible conclusion that while the long period and the regular 



behaviour of the modulation of RZ Lyr favour the occurrence 
of modulation-phase locked realization of the resonance con- 
dition, in RR Lyr, which has a shorter Blazhko period and 
its consecutive modulation cycles are not fully repetitive, 
the connection between the appearence of period doubling 
and modulation phase is only partial if there exists any. 

• The modulation of RZ Lyr is intrinsically doubly peri- 
odic. The secondary modulation period is close to but does 
not exactly equal with one quarter of the period of the 121-d 
main modulation. 

Many Blazhko stars exhibit complex, not strictly periodic 
light-curve modulation. The modulations of these variables 
are most probably multi-periodic. The only detailed study of 
a double-modulation-period Blazhko star, CZ Lac, showed 
that there is a tendency that the periods/frequencies of 
the independent modulat ion comp onent s are close to be- 
ing in resonance (|S6dor et al.ll2011l ). This is also the case in 
RZ Lyr, where / m i//m2 w 1 : 4. It is also worth to note that 
the two modulations have different characters both in RZ 
Lyr and in CZ Lac. The explanation of the multi-periodic 
modulations is one of the greatest challenge for the theoret- 
ical elucidation of the Blazhko effect. 

• By the aid of the IP method, the changes in the mean 
values of the physical parameters during the modulation are 
determined. An extremely large, ~3 L© luminosity variation 
is detected; the star is the brightest at the large-amplitude 
phase of the modulation. The variations of the physical pa- 
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rameters connected to the primary and the secondary mod- 
ulations show different behaviours. The pulsation ampli- 
tude, period, stellar radius, luminosity and temperature all 
vary in phase with each other during the primary modu- 
lation. On the contrary, during the secondary modulation, 
the pulsation amplitude and period variations have oppo- 
site directions, and the radius variation follows the period, 
while the marginal luminosity variation follows the ampli- 
tude changes. 



According to the idea of lStothersI (|200rj . l201Ct ) a convective 
dynamo cycle is responsible for the light-curve modulation of 
Blazhko stars. The model gives predictions for the strengths 
of the period variation during the modulation as well as for 
the phase relation of the period and radius variations. No 
phase (period) variation of variables with pulsation period 
around 0.5 d, and anticorrelated and correlated changes in 
the pulsation period and the radius variations are predicted 
to occur in hotter and cooler variables, respectively. The ob - 
served physical variations of MW Ly r (|jurcsik et al.|[2008bl ) 
and DM Cyg (|jurcsik et al . 2009a) agree well with these 
pred ictions, while the complex modulation behaviour of CZ 
Lac (ISodor et al 1120111 ) fit in only partially with this scenario 
|Stothersll201ll 



The 0.51-d pulsation period of RZ Lyr is close to the 
half-day critical period value, and indeed, its amplitude vari- 
ation is much larger than the phase variation during the 
primary 121-d modulation. It is, however, not true for the 
secondary modulation, which is dominated by phase varia- 
tion. Note that there are Blazhko stars showing negligible 
phase variations in the 0.55-0.60 d pulsation period regime, 
too (Konkoly Blazhko Survey II, Jurcsik et al. in prepara- 
tion). 

The phase relations between the period and radius 
changes are correlated with (dP/P)/(dR/R) ratios of 
0.0011/0.006 and 0.0018/0.002 for the primary and sec- 
ondary modulations, respectively. RZ Lyr is found to be 
indeed hotter than MW Lyr, which shows anti correlated pe- 
riod and radius changes \ Jurcsik et al.ll2008bl ). but its tem- 
perature is the same or e ven somewhat warm er than the 
temperature of DM Cyg [1 Jurcsik et al.l l2~009al ). a smaller 
amplitude 'twin' of MW Lyr. 

These results confirm only in p art Stother's mode l, 
which has been seriously criticized bv lSmolec et~aH (|201ll ). 
The large variety of the physical parameters of RR Lyrae 
stars makes it, however, unlikely that any model could give 
universally valid, unique predictions that hold for the whole 
possible metallicity and luminosity ranges of the variables. 
Mapping the physical parameter variations of a larger sam- 
ple of Blazhko stars during the modulation cycle is required 
to reveal any connection between the amplitude and phase 
relations of these changes and the physical parameters of the 
stars. 
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Table 9. Konkoly CCD observations of the two new variables 
discovered in the field of RZ Lyr. The full table is available as 
Supporting Information with the online version of this paper. 
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Figure 13. BVIq and B — V, V — Iq light and colour curves of 
VI, a fi Lyr-type new variable in the field of RZ Lyr. 

8 APPENDIX 

As given in Section 12.11 two new variables in the field of 
view of RZ Lyr have been identified. The photometric data 
of the new variables are given in the electronic edition of 
this paper (Table [5J). VI and V2 are measured relative to 
the comparison stars 'a' and 'c', respectively (Fig. [TJ. The 
light and colour curves these variables are shown in Figs. [13] 
and 1141 respectively. Basic information on VI and V2 are 
summarized in Table [8] 
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Table 8. Basic information on the new variables 



Variable 


2MASS ID 


(V) 


(B-V) 


(V - Ic) 


Type 


Period 


T [HJD] 


VI 


18431843+3248587 


11.82 


0.23 


0.28 


EB 


0.927033(6) 


2455353.390 


Y2 


18435377+3247459 


14.19 


0.77 


0.86 


EA 


1.447245(6) 


2455396.370 




0.5 



phase 



Figure 14. BVIq and B — V , V — Iq light and colour curves of 
V2, an Algol-type new variable in the field of RZ Lyr. 



